About Landau-Hopf scenario in a system of coupled self-oscillators 
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The conditions are discussed for which an ensemble of interacting oscillators may demonstrate 
the Landau-Hopf scenario of successive birth of multi-frequency quasi-periodic motions. A model is 
proposed that is a network of five globally coupled oscillators characterized by controlled degree of 
activation of individual oscillators. Illustrations are given for successive birth of tori of increasing 
dimension via quasi-periodic Hopf bifurcations. 
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I. INTRODUCTION 

In this paper we discuss a problem that attracts at- 
tention for many years, but has not received yet a sat- 
isfactory solution and understanding. It consists in a 
question of possibility of observation of the scenario of 
transition to complex dynamical behavior suggested by 
Landau and Hopf, e.g. in the context of hydrodynamic 
turbulence. This scenario includes successive secondary 
quasi-periodic Hopf bifurcation, which follow the primary 
one, being accompanied with a sequential birth of incom- 
mensurable frequency components of the motion. We 
suggest here some novel approach to this long-standing 
problem based on artificial constructing and examina- 
tion of model systems with the required type of behavior, 
which allow physical implementation (say, in mechanics 
or electronics). A concrete example we propose is com- 
posed as a network of coupled van der Pol oscillators. 
It may be thought as a setup universal in some physi- 
cal sense, due to the fact that the van der Pol oscillator 
is associated with a normal form for the Andronov-Hopf 
bifurcation. 

Back in the forties of the last century, Landau [l| , and 
later Hopf 0, suggested that in the course of variation of 
the control parameter (Reynolds number) the hydrody- 
namic turbulence may appear via the following sequence 
of events. First, after the stability loss of a steady flow, a 
self-oscillatory regime appears that corresponds to a limit 
cycle in the phase space of the system. (This transition 
is called now the Andronov-Hopf bifurcation.) Then, in 
turn, the periodic motion becomes unstable and under- 
goes the secondary bifurcation. (In modern terminology 
it is referred to as the Neimark-Sacker bifurcation.) Now, 
the dynamics takes place on a two-dimensional torus in 
the phase space of the system, born as a result of the 
soft bifurcation transition. Next, this motion under- 
goes a similar bifurcation, resulting in emergence of one 
more oscillatory component, with its own independent 
frequency, giving rise to a three- frequency torus, and so 
on. As a result of successfully growing number of the 
incommensurable frequencies in the spectrum of oscilla- 
tions, the regime is complicating, and finally the turbu- 
lence develops. 

An attractive feature of the Landau argumentation is 



its generality; in fact it does not appeal to a specific 
physical nature or to concrete mathematical equations of 
the system, and, hence, must have equal relation to a 
variety of multidimensional and spatially extended dissi- 
pative nonlinear systems. 

Later the approach of Landau-Hopf was criticized in a 
frame of the concept of Ruelle and Takens [3|] , who indi- 
cated a possibility of destruction of low-dimensional tori 
due to arbitrarily small variations in the system evolu- 
tion operator accompanied with the onset of a strange 
attractor corresponding to chaotic dynamics. 

To verify the Ruelle- Takens concept, several re- 
searchers performed numerical experiments for model 
systems and observed that the low-dimensional tori usu- 
ally survive under small perturbations, although these 
tori may be destroyed givingbirth to chaos at sufficiently 
large magnitude of them @,[5|. Hence, a widespread con- 
clusion of the universality of the Ruelle- Takens scenario 
should be regarded as hasty and ill-founded. 

On the other hand, one can criticize the Landau-Hopf 
picture in a frame of the concept of synchronization, 
which is a universal nonlinear phenomenon [6]. Indeed, 
in situations of close basic frequencies for the involved 
oscillatory components (or of some combinations of their 
frequencies), the system may demonstrate either attrac- 
tive periodic orbits, or lower dimension attractive reso- 
nant tori placed on higher-dimensional tori [7MTll| . So, 
resonances and the possible occurrence of chaos are ex- 
pected to modify essentially the picture of dynamical 
regimes and bifurcations; it makes the Landau-Hopf sce- 
nario problematic. In this context it would be interesting 
nevertheless to search for specific systems, in which such 
a scenario actually could occur, at least, one could ob- 
serve a sufficiently large number of initial quasi-periodic 
bifurcations. Constructing an appropriate model based 
on a set of coupled nonlinear oscillators is the purpose of 
the present article. 

When choosing a model to analyze, we take into ac- 
count several important physical aspects. First, the os- 
cillatory modes responsible for the dynamics must be 
characterized by different degrees of activation. If we 
consider, say, an ensemble of coupled van der Pol oscilla- 
tors, we must introduce a set of parameters controlling 
the negative friction and select them appropriately to en- 
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sure conditions for gradual involving of the modes in the 
motion in the course of decrease of dissipation level. (In 
the context of the hydrodynamic problems it just corre- 
sponds to increasing Reynolds number.) Second, all the 
oscillatory modes have to be separated in frequency in 
sufficient degree. Otherwise, essential interaction of the 
modes would occur, which can destroy the Landau-Hopf 
picture. Third and finally, it is desirable to have a situ- 
ation, where the coupled oscillators are arranged in such 
way that they all are involved essentially in the interac- 
tion; it means that no preferable interaction of each con- 
crete partial oscillator should occur, say, with spatially 
close neighbors, or with different number of the neigh- 
bors. In such situation the control parameter of each 
single oscillator A^ will regulate a certain quasi-periodic 
bifurcation. Otherwise, it might happen that due to the 
different number of interacting neighbors the elements 
with close level of activation would have essentially dif- 
ferent influence on the picture of the multi-frequency dy- 
namics fn\ . 

Selecting the model, it is crucial to guarantee the oc- 
currence of the definite type of the bifurcations. For the 
Landau-Hopf scenario these must be bifurcations of soft 
(non-catastrophic) birth of tori of increasing dimension. 
(In modern literature such bifurcations are usually re- 
ferred to as the quasi-periodic Hopf bifurcations [~L2lfl3|.) 
In contrast, resonance mode locking phenomena are as- 
sociated with saddle-node bifurcations of tori. Unfortu- 
nately, no simple numerical algorithms for identifying the 
quasi-periodic bifurcations are elaborated. Qualitatively, 
the types of bifurcations of quasi-periodic regimes can be 
identified in computations by means of analysis of the 
behavior of the Lyapunov exponents, and this method 
will be used henceforth. 



II. NETWORK OF FIVE COUPLED 
SELF-OSCILLATORS 

Let us construct a model satisfying the above require- 
ments. It will be a kind of network of five oscillators 
with global coupling (that implies coupling of each el- 
ement with each other). We suppose the oscillators to 
have equidistant spectrum of natural frequencies and 
leave a single frequency parameter A in the model, which 
controls the mutual mismatch between the oscillators. 
For large values of this parameter the oscillators will be 
desynchronized. We assume that all oscillators are char- 
acterized by different values of the parameters respon- 
sible for exceeding the excitation threshold; it will al- 
low gradual activation of the corresponding modes in the 
course of the transition under consideration. 

The set of equations for the model system we propose 
reads 



x — (Ai — x 2 )x + x + j(4:X — y — z — w — v) = 0, 

y-{\2- y 2 )y + (l + j)y + £(4tfr -±-z-w-v) = o, 

z-(X 3 - z 2 )z + (1 + f)z + f (4i -y-x-w-v)=0, 
w - (A 4 - w 2 )w + (1 + + %(4w -y-x-z-v) = Q, 

v - (A 5 - v 2 )v + (1 + A)v + f (4v - ij - x - z - w) = 0. 

(1) 

Here A^ are control parameters responsible for exci- 
tation of the partial oscillators, A determines the fre- 
quency detuning of the oscillators, and the frequency 
of the first oscillator is unit. We set hereafter Ai = 
0.1, A 2 = 0.2, A 3 = 0.3, A 4 = 0.4, A 5 = 0.5. 

Our main illustration will be a Lyapunov chart on the 
parameter plane for frequency detuning A and dissipative 
coupling \i. To draw this chart we proceed as follows [l(J 
1 1 II ] , First, select a point on the parameter plane and 
compute the spectrum of Lyapunov exponents for the 
Poincare map of the system (pQ) there. Then, analyzing 
the largest four exponents, take into account their sings 
and evaluate a number of zero exponents among them to 
determine a type of the dynamical regime (the attractor 
type) in the system. Namely, it is determined as 

1. a limit cycle P, if Ai < 0, A 2 < 0, A 3 < 0, A 4 < 0, 

2. a two- frequency torus T 2 , if Ai = 0, A 2 < 0, A3 < 
0, A 4 < 0, 

3. a three- frequency torus T3, if Ai = 0, A 2 = 
0, A 3 < 0,A 4 < 0, 

4. a four- frequency torus T 4 , if Ai = 0, A 2 = 0, A3 = 
0, A 4 < 0, 

5. a five- frequency torus T5, if Ai = 0, A 2 = 0, A3 = 
0, A 4 = 0, 

6. chaos C, if Ai > 0, A 2 < 0, A 3 < 0, A 4 < 0. 

Then, the respective pixel on the parameter plane is 
attributed with a definite color, in accordance with the 
detected dynamical regime. Scanning the entire parame- 
ter plane, we get the two-parameter chart of the dynam- 
ical regimes. 

Figure Q] shows the Lyapunov chart obtained for the 
model (pQ) on the parameter plane (A, ji). Because of the 
structure of the system, the main resonance effects are ac- 
tually excluded. In particular, no tongues of resonant tori 
of different dimensions similar to those mentioned e.g. in 
Refs. El are observed. The only pronounced tongue 
corresponding to a resonant two-frequency torus may be 
seen in the region A < 0.5, between the domains of the 
five- frequency tori and of the complete synchronization. 
Interestingly, it is immersed in the region of chaos that 
occurs at small coupling [16|. Here the three- frequency, 
four-frequency, and partially the five- frequency tori are 
destroyed, although the chaos is weak. (Typical values of 
the positive Lyapunov exponent are of order lO^-lO -2 
here.) To some extent, one can say that the Ruelle- 
Takens scenario occurs in this region. Note additionally 
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FIG. 1: Chart of Lyapunov exponents for the network of five 
globally coupled identical van der Pol oscillators (pQ), where 
Ai = 0.1, A 2 = 0.2, A 3 = 0.3, A 4 = 0.4, A 5 = 0.5. 

that a small island of chaos takes place surrounded by 
the area of two-frequency tori, at the coupling coefficient 
values about 0.3-0.4. Hence, to give rise to chaos it is 
not needed to deal necessarily with destruction of three- 
frequency tori; it appears to be sufficient to have the 
smaller dimension. 

On the other hand, in the case of large detuning of the 
oscillators, A > 1, with decreasing dissipation parame- 
ter jjL one can observe the appearance of all the tori of 
higher dimension. Boundaries of the relevant areas in 
the asymptotic A — > oo correspond to the values of the 
control parameters /i=A^. (For convenience, these values 
are indicated by arrows on the right edge of the chart.) 
Thus, a decrease in the parameter of dissipative coupling 
qualitatively corresponds to the pattern expected for the 
Landau-Hopf scenario. 

III. THE CASCADE OF QUASI-PERIODIC 
BIFURCATIONS 

To be confident that we deal here really with a cas- 
cade of quasi-periodic Hopf-Neimark-Sacker bifurcations 
for tori of higher dimensions, we turn to the graphs of the 
Lyapunov exponents for the Poincare map of the system 
in Fig. [2j The exponents are plotted versus the cou- 
pling parameter /i, along the vertical line A = 3 in the 
parameter plane. One can see that for large values of 
the coupling parameter, all the exponents are negative 
that corresponds to the so-called oscillator death regime 
(OD). This mode is typical for coupled oscillators, and 
occurs due to the dissipative nature of coupling: at suf- 
ficiently strong couplings it dampens the oscillations 
The point H corresponds to the Andronov-Hopf bifur- 
cation, and here one exponent Ai becomes equal zero. 
Then, at the Neimark-S acker bifurcation NS, one more 
exponent A2 vanishes. At the point QH\ the next ex- 




FIG. 2: Graphs of six Lyapunov exponents of the Poincare 
map of model (pQ) at the value of the detuning parameter A=3. 

ponent A3 becomes zero that corresponds to the birth 
of the three-dimensional torus. To identify the nature of 
the bifurcation, we note that before the bifurcation the 
exponents A3 and A4 are equal. This is a characteristic 
feature of the soft quasi-periodic Neimark-S acker bifur- 
cation [12|. Alternatively, at the quasi-periodic saddle- 
node bifurcation the Lyapunov exponents behave differ- 
ently m d3 . 

Then, exactly in the same manner the remaining expo- 
nents behave: at the point QH2 a four-dimensional torus 
appears due to the secondary quasi-periodic Neimark- 
Sacker bifurcation, and at the point QH3 a five- 
dimensional torus arises. 

Figures [3] and 3] illustrate the development of the 
Landau-Hopf scenario. 

Fig. [3] demonstrates evolution of the phase portraits 
of the oscillators with decrease of the dissipation pa- 
rameter. Each diagram shows portraits for all five os- 
cillators on their phase planes (the generalized coordi- 
nate versus the generalized velocity), and the orbits re- 
lating to the different oscillators are drawn in different 
colors. Fig. [3^ corresponds to a large enough value of 
the dissipation parameter /i = 0.45. In this case we have 
a limit cycle embedded in the multi-dimensional state 
space of the system. It is seen clearly that the oscil- 
lations of the fifth oscillator are the most intense (it is 
so because the control parameter A of this oscillator is 
the largest). The intensity of the motion for the other 
oscillators decreases. In panel (b) corresponding to the 
three- frequency torus at \i = 0.25, the trajectories of all 
oscillators are perturbed quasi-periodically. However, in 
panel (c), at \i = 0.05, where the five- frequency torus oc- 
curs, the magnitude of the perturbations falls again: the 
small coupling weakly perturbs the orbits. 

Figure |4] shows evolution of the Fourier spectrum of 
the fifth oscillator in the course of development of the 
Landau-Hopf scenario in the system (pQ). One can ob- 
serve the gradual enrichment of the spectrum with new 
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FIG. 3: Evolution of the phase portraits of oscillators of the 
system (pQ) with decrease of the dissipation parameter: a) 
[i — 0.45, the limit cycle, b) /x = 0.25, the three-frequency 
torus, c) ijl , = 0.05, the five- frequency torus. The numbers 
numerate partial oscillators of the system; the frequency de- 
tuning parameter is A = 3. 



spectral lines corresponding to the quasi-periodic regimes 
of the increasing dimension. Note meanwhile that in the 
panel (e) the degree of complexity of the spectrum visu- 
ally rather decreases, although, strictly speaking, the di- 
mension of the motion grows (the five-dimensional torus). 
The reason is the significant reduction in the level of cou- 
pling: at very low coupling the attractor is perturbed 
very weakly, and is close to the limit cycle of the au- 
tonomous oscillators (panel (c)). Therefore, the height 
of the corresponding spectral lines decreases, although in 
smaller scales the spectrum remains complex. 

Figure [5] shows a plot for the observed frequencies ver- 
sus the coupling parameter for the case A = 3. The fre- 
quencies are defined as the average rates of variation in 
time for the phases of the partial oscillators: uji =< (pi >. 
They are evaluated numerically by counting a number 
of crossings of a certain surface in the phase space for 
a long time interval. The surfaces selected correspond 
to zero values for coordinate variables for the respec- 
tive oscillators. From Figure [5] one can see that in the 
range 0.4 < \i < 0.5, for the coupling level large enough, 
all the frequencies are the same, approximately equal to 
the frequency of the most excited fifth oscillator. With 
the decrease of the coupling strength, at fi « 0.4, a 
Neimark-S acker bifurcation occurs; in Figure [5] one can 
observe the appearance of one more frequency, which is 
close to the natural frequency of the fourth oscillator. 
Next, at \i « 0.3, one more new frequency branches out, 
which quickly approaches a value corresponding to the 
natural frequency of the third oscillator. A specific in- 
teresting behavior occurs for the remaining two oscilla- 



tors. Their frequencies undergo a sharp change almost 
simultaneously with the third oscillator, but in the range 
0.2 < [i < 0.3 they follow approximately the higher fre- 
quency of the fourth oscillator. In other words, the ex- 
cited mode forces the weakly excited oscillators to adjust 
their frequencies to the former one. Then, the fourth 
frequency branches out, and, finally, at the parameter 
value ji « 0.1, the fifth frequency appears. [17| As we see, 
the picture of the "tree of synchronization" arises in the 
system manifesting the sequence of quasi-periodic Hopf- 
Neimark-S acker bifurcations. This tree is of a rather spe- 
cific sort, essentially distinct from that intrinsic to the 
phase oscillatory systems of traditional kind (see, for ex- 
ample, Figure 1.11a in [6], or Figure 1 in [15], etc.). In 
our case, one of the frequencies remains almost constant 
in the course of each bifurcation, while a new component 
branches out sharply. This is so because of the differ- 
ence is the involved types of bifurcations: in Refs. [H, [TBJ 
these are the saddle-node bifurcations of invariant tori, 
while in our case these are the secondary quasi-periodic 
Hopf-Neimark-Sacker bifurcations. 



A. Conclusion 

Thus, under certain conditions, e.g. in the case of non- 
identical parameters of the active modes and of detuning 
for the modes in ensemble of self-oscillating elements, one 
can observe a sequential cascade of soft quasi-periodic bi- 
furcations involving tori of increasing dimension that can 
be regarded as occurrence of the Landau-Hopf scenario. 
The model outlined in this paper can be implemented, 
say, as an electronic device, which will serve as an exam- 
ple of a physical system that demonstrates this scenario 
under variation of the control parameter. 
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FIG. 4: The evolution of the Fourier spectrum fifth oscillator 
in (pQ) when the dissipation in the system, A = 3, a) /i = 0.45, 
b) fj, = 0.35, c) ii = 0.25, d) ii = 0.15, e) ji = 0.05. 



